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ABSTRACT 

We investigate the Dijkgraaf-Vafa proposal when supersymmetry is broken. We consider 
U(N) SYM with chiral adjoint matter where the coupling constants in the tree-level su- 
perpotential are promoted to chiral spurions. The holomorphic part of the low-energy 
glueball superpotential can still be analyzed. We compute the holomorphic supersymme- 
try breaking contributions using methods of the geometry underlying the M = 1 effective 
gauge theory viewed as a Whitham system. We also study the change in the effective 
glueball superpotential using perturbative supergraph techniques in the presence of spu- 
rions. 
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1 Introduction 



Dijkgraaf and Vafa have proposed that the low-energy glueball effective superpotential 
of M = 1 supersymmetric gauge theories in four dimensions can be computed via an 
auxiliary matrix model [T]. The simplest case is a U(N) gauge theory coupled to a 
massive adjoint chiral matter multiplet $ with a tree-level superpotential W(&). The 
proposal stems from a set of string dualities in the framework of geometrically engineered 
gauge theories, topological strings and matrix models (21 El E] ■ The large-N matrix model 
analysis brings in an algebraic curve which may correspond to a Calabi-Yau dual geometry 
[2]. We shall consider gauge theories that can be obtained from string theories that lead 
to such geometries. The DV proposal has been tested and supported directly on the field 
theoretical side by perturbative computation via superfields formalism [3] and then by 
using arguments based on anomaly equations [Hj. 

We study here the case where susy is broken explicitly (soft and/or non soft) by 
the introduction of spurionic fields [(i . Holomorphy at large is lost, but holomorphic 
quantities such as the glueball superpotential can be still analyzed and one can compare 
the computation in the superfields formalism adapted to spurion fields with that one 
using the algebraic curve underlying the effective gauge theory. In order to discuss such 
breaking we utilize two notions: a closed string realization of the method of the spurions 
[7j and Whitham deformations [HI E] . 

It is worth to recall the geometrical origin of such gauge theories for type IIB string 
theory in order to insert the notion of spurion in a natural way in this language. We have 
in mind D-branes partially wrapped over non trivial 2-cycles of non compact CY and the 
dual description where D-branes have been replaced by fluxes [2j. In the UV, adjoint 
chiral multiplets $ arise from holomorphic deformations of the supersymmetric cycles 
and of open string gauge bundles on these cycles. A four- dimensional superpotential for 
these fields can arise and can be written as W = W(<&,gk) where depend only on 
the complex structure. From the perspective of the .D3-brane action in the low-energy 
limit, where supergravity decouples, the g% can be interpreted as couplings. As already 
suggested in [7j, the susy breaking parameters are described by auxiliary components of 
the closed string fields, typically magnetic fluxes along CY directions, depending on the 
complex structure moduli. Such fluxes are introduced by hand without back reaction 
of the string or of the supergravity backgrounds. In the four dimensional supergravity 
language they are F-components of chiral multiplets which depend only on the complex 
structure moduli. Vev of such F-terms cause spontaneous breaking of local susy and, in 
the appropriate flat limit with decoupling of supergravity, they appear as explicit breaking 
terms which can be written in the spurionic fashion in the rigid susy action. 

A non-perturbative analysis of susy broken effective dynamics has been done in ^U] 
for M = 2 supersymmetric gauge theories. In that context the connection between the 
Seiberg-Witten solution jTT] and integrable systems (Whitham hierarchy) |H. was used. 
The authors of JUj break susy promoting the Whitham parameters of the hierarchy to 
spurions and then compute the broken effective potential using the M = 2 integrable 
structure. 
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As in the M = 2 relation between the Whitham systems and the TV = 1 effective 

geometry was established in jS]. This suggests to break supersymmetry promoting the 
Whitham parameters to spurions as in the M = 2 case. In the M = 1 geometry the 
Whitham parameters are precisely the tree-level coupling costants of the matter super- 
potential 9\. We will break the M = 1 supersymmetry promoting them to spurions, and 
the Whitham hierarchy can then be interpreted as a family of supersymmetry breaking 
deformations of the original theory. Using this interpretation, we will compute directly 
from the geometrical data the holomorphic supersymmetry breaking contributions in the 
low-energy effective glueball superpotential. 

We have also analyzed with perturbative supergraph techniques the effective glueball 
superpotential when susy is broken with spurions. Arguments for the computability of 
the effective superpotential have been presented in [7j. If supersymmetry is broken, holo- 
morphicity in the coupling constants is no longer guaranteed, the computation is much 
harder than in the M = 1 case and the simplifications of || do not work in general. 
Anyway, we can restrict ourselves to a particular subclass of contributions for which a 
spurionic superfields generalization of the techniques in jl] can be done. Within such 
strong approximation and with unbroken U(N) gauge group, we find that to all order in 
the glueball superfield the effective superpotential has the same functional form of the 
J\f = 1 case where the coupling constants are replaced by spurions and so it results still 
holomorphic. 

The paper is organized as follows: In section 2 we review the geometry underlying the 
Dijkgraaf-Vafa proposal. In section 3 we introduce supersymmetry breaking by spurions 
and discuss the low-energy glueball superpotential. In section 4 we discuss the geometry 
as a Whitham system and use it in the susy broken case. In section 5 we treat the explicit 
example of a deformed susy broken cubic tree-level superpotential. In section 6 we use 
perturbative superspace techniques along the line mentioned above. Section 7 is devoted 
to conclusions. At the end, in two appendices, we describe the computational details of 
section 5 and section 6. 

2 The geometrical picture 

We consider the particular case of a Af = 1, U(N) gauge theory with a degree n + 1 
polynomial tree-level superpotential for the chiral matter superfields in the adjoint 

representation of the gauge group 



In a generic vacuum the gauge group U(N) is broken to U(Ni) x • ■ ■ x U(N n ). In the IR 
limit the effective low-energy degrees of freedom are described by the glueball superfields 
Si = -^sTrW^Wai where is the fermionic chiral superfield, field strength of the 
vector multiplet of the unbroken gauge group U(Ni). 



n+1 




(1) 



k=l 
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The expression for the non perturbative glueball superpotential reads 



w.//(s) = -E 



NiT^ + 27iiTiSi 



(2) 



where J 7 is the prepotential which can be computed from the geometrical data [3 E] • I n 
[l] it has been proposed to reinterpret and compute this prepotential as the free energy of 
an associated matrix model. In H2] it was also deduced directly on the field theoretical 
ground using generalized Konishi anomaly equations. 

The geometry associated with the low-energy theory is described by a family of genus 
g = n — 1 Riemann surfaces and by a meromorphic differential dS 

V 2 = [W(x)f + f^~ 1 \x), (3) 
dS = ydx= \/[W(x)f + f( n -V(x)dx . (4) 

The degree n — 1 polynomial f^ n ~ l \x) = Ym=o fi xl 5 i s associated with the quantum 
contributions and the coefficients fi (I = 0, • ■ • , n — 2) are the moduli of the complex 
curve; the derivatives of the meromorphic differential (jlj with respect to the moduli gives 
holomorphic differentials. 

A basis of canonical cycles |HJ H2] is {a\ Pi, a , /3 }, where i — 2, . . . , n, with intersec- 
tion numbers (/3b Ha" = 5%). The cycles are all compact except (3q. We label the cuts 
starting from the larger real root of the algebraic curve (jBJ), so from right to left. The 
a*-cycle surrounds counterclockwise the i-th cut while the a°-cycle encircles all the cuts 
and then gives the residue at infinity. The dual /3j-cycle {i = 2, . . . , n) passes clockwise 
through the i-th and the first cut, while flo goes from the second sheet infinity to the 
first passing through the first cut. The periods Si, the parameter t and the conjugated 
periods are defined as 



Si = (b dS 



t = <f dS = -Res^dS) = / n 1 , (5) 
YU = \<f> dS , n = ~ / dS . (6) 



or — ' u 9 , 



In these variables the effective superpotential computed by the geometry is 



n dJ- n dJ- 

-W eff = NU + J2N l U t = N— + J2N l g- , (7) 

i=2 i=2 * 



where Ylj=i Nj = N. In the previous formula we have introduced the prepotential 4 JF 
such that its derivatives w.r.t. the {si,t } periods give the dual ones {IL,I1 }. 

4 The prepotential differs from the usual one ^ |2] for a multiplicative factor due to the change of 
variables. Anyway, this difference is not felt by the effective superpotential We// which is a function of 
the dual periods, the quantities which really enter in the computation, as in i|7|l. 
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Upon getting the superpotential as a function of the variables and to, we return to 
the variables of 12] using 5 

Si = 2>Sj , i = 2, . . . , n , 

n 

to = -2^2 Si , (8) 
in fact the Si are the physical variables which are interpreted as the glueball superfields. 

3 Super symmetry breaking 

The introduction of spurionic fields provides the standard mechanism for the explicit (soft 
and/or non soft) breaking of global supersymmetry. In the M = 1 case the tree-level 
superpotential W tree , and the effective glueball prepotential J 7 , depend on the coupling 
constants g m associated with the operators Tr $ m in the ultraviolet action. In order 
to break M = 1 supersymmetry down to M = we promote the coupling constants 
g m to M = 1 chiral superfields G m and then we freeze the scalar and the auxiliary F- 
components to constant values. In this way the chiral spurions G m = g m + 9 2 F m produce 
non supersymmetric terms in the superpotential W tre e- We want to study their effects 
on the low energy glueball effective superpotential under the assumption that the low 
energy degrees of freedom are still the glueballs. The breaking parameters T m must be 
considered the smallest scales in the theory. They are thought as small perturbations of 
the M = 1 theory by keeping fixed the M = 1 vacuum structure and the gauge symmetry 
breaking patterns U(N) -> U(Nx) x • • • x U(N n ). 

We set the scalar components of G m equal to the coupling constants g m for m < n + 1, 
zero for m > n + 1, and the F-components T m will be considered as small susy breaking 
parameters for all G m . Explicitly 

G k = g k + 9 2 T k , k<n + l , (9) 
Gj = e 2 Tj , j > n + 1 , (10) 

and hence we will consider tree-level superpotential ((IJ perturbed as 

n+l „ „ 

w tree (<5>) = Y.T Tl ® k + e2 S Tr$ ' • ( u ) 

fc=l j>n+l 3 

Notice that besides having promoted to spurion the coupling constants already appearing 
in the tree-level superpotential, we have also added pure auxiliary F-terms. For k > 3 
these spurionic terms are not soft and quadratic divergences can appear in the wave 
function renormalization; in any case they have to be considered as dangerously irrelevant 

5 The variables of are Sj — — \ § Aj dS and IT, = \ § B . dS, with {AjjB^; j = 1, . . . ,n} a different 
set of cycles with all Bj non compact. 
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operators with the usual warning [T3*l [7j. The T m for m < n + 1 can be interpreted as 
vacuum expectations values of fluxes [7j, whereas it is not obvious that this is the case 
for m > n + 1. In any case, we will see that the generalization to all the T m terms is of 
some interest in the application of the Witham approach. 

Let now analize what happens in the effective theory when the T m are turned on. 
We will restrict ourselves to a discussion of some formal aspects which can be extracted 
from the geometry of the M = 1 case. We assume that in the effective dynamics the 
emergence of the spurions G m are controlled by the holomorphic dependence of the M = 1 
prepotential F(Si,g m ) on the coupling constants. If we restrict ourselves to holomorphic 
terms in the low-energy glueball superpotential, the prepotential in the susy broken phase 
has the same functional form as the M = 1 case where now the coupling constants g m 
are replaced by the spurions as G m . This is essentially a naturalness assumption on 
the effective superpotential [Hj. In section 6 we will discuss these assumptions using 
superfields perturbative techniques extending jl] to the susy broken case. 
We make some comments about the interpretation of the couplings Tj (j > n + 1). They 
must be understood as coming from tree-level superpotential W tre e of degree greater than 
n + 1 where also the scalar coupling constants gj above the (n + l)-degree are turned on. 
The low energy glueball prepotential will also depend on all these couplings. We then 
consider the effective theory of (jllj) as obtained from that one of higher degree in the limit 
where gj — > (j > n + 1) and in the same vacuum of the theory of (n + l)-degree 6 . In 
conclusion the prepotential depends on the n glueball superfields Si (in our conventions 
t and Si) and it is evaluated where gj =0. 

We expand now the prepotential J-{S i: G m ) around the supersymmetric vacuum. If 
we consider the case of broken supersymmetry with G m having the form (0 HO]) the terms 
with more than one power of r m will not give any contribution and we have 

^(s^g^T^Tj) = F(si, g m )\ g .= + 



The first term in this expression is the prepotential of the supersymmetric case for a 
theory with tree-level superpotential of (n + l)-degree. As just discussed the last term 
is interpreted as coming from an higher degree theory in the appropriate limit. 
We now insert this expression in (jjj) and we obtain the holomorphic glueball superpotential 
associated with a tree-level susy breaking superpotential as (|TT|) 




(12) 




(13) 



6 We can choose J\f = 1 massive theories with classical vacua configuration which are nonsingular for 
gj — > such that there is analitycity of the glueball superpotential around gj = 0. 
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In this expression the first terms within the square bracket are supersymmetric, whereas 
the others break susy explicitily: they involve second derivatives of the prepotential eval- 
uated where the gj = (j > n + 1). We will show in the next section how to obtain 
directly and efficiently from the geometrical data of the M = 1 theory the mixed sec- 
ond derivatives of T appearing in (fT3|) in order to extract the effective supersymmetry 
breaking contributions. 



4 The J\f = 1 geometry and Whitham systems 

The geometry of the M = 1 low-energy effective theory is associated with the generating 
meromorphic differential dS and it can be thought as coming from a Seiberg-Witten 
geometry of a Af = 2 theory [3]. The addition of a superpotential together with a 
geometric transition and a desingularization leads to such geometry with parameters g k 
and complex moduli f] 2 . The couplings g k can be viewed as Whitham deformations 
of the previous SW geometry. Performing a Whitham deformation mean extending the 
parameter space of the curve with extra variables [S] . As a consequence of this deformation 
the moduli of the curve and also the generating differential become functions of these new 
parameters. 

As shown in jHj the M = 1 geometry can be embedded into the Whitham framework. 
The moduli fi of the curve © are functions fi = fi(gk, to, s») of the Whitham parameters 
g k and of (£ , the periods of the generating differential dS along the a-cycles. We 
review some results of and set up our conventions. 

One of the advantages we gain using Whitham description is that it provides an effi- 
cient way to compute the mixed second derivatives appearing in (|13|) directly in terms of 
geometrical data since the coupling constants are considered as independent parameters. 

Using the whole set of variables (gk, t , si) characterizing the curve Q and the gener- 
ating differential (J1J), the Whitham system can be defined by the following set of equations 

— — = duji , — — = dllo , — — = dllk , (14) 
dsi dt dg k 

where duji are normalized holomorphic differentials 

ddS 



dujj = 5ij . (15) 



The differentials dQk are meromorphic of the second kind with poles only at the infinity 
points ±oo; dQ is a differential of the third kind with residue at ±oo. They have vanishing 
a-periods and behave at infinity as (£ = -) 

/ dn = |^ = o , / dn k = p- = o ; tzn, = -(c 1 ' 1 + o(i)K • (16) 

J a i dt J ai dg k 

These normalization conditions characterize Sj, to and g k as independent variables. 
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The generating differential dS is then a linear combination of the differentials ()14j) 



n+l 



dS = Si duJi + t dClo + g k dfl^ 

i=2 fe=l 



\ 



n-2 



[Wixy^ + J^h^ + ^n+itox^dx . (17) 



fe=0 



Consinstency of the equality in (|T7j) requires that 

g k = -Re Soo+ (x- k dS) , (18) 
which can be verified j^j. Using (J 17)) the meromorphic differentials dQi can be written as 



dn - ddS - 9n+l%n ~ l dx i l ^ dfixl dx 

dQ, - k-l n 

ogk V 2 ^ dg k y 

ddS [W'(x)x n + t x n - 1 ] J 1^ 9// x 1 

cft2 n+1 = = dx +n/,a dx ■ ( 19 ) 

og n +i y 2 ^ ^ n+ i y 

In this framework, the prepotential T and so the special geometry can be introduced 
thanks to the Riemann bilinear relations which guarantee the integrability condition of 
the prepotential J2|. We must define correctly the first derivatives of T with respect to 
both the periods and the coupling constants 

dT „ 1 f ,„ dT „ 1 r°° + ,„ dT „ (x k 



8si ^2 t dS ■ wr n °-2L dS - wr^W 3 ) ■ (20) 

As we have seen in the previous section, the supersymmetry breaking contributions ap- 
pearing in the effective glueball superpotential (fTHj) are mixed second derivatives of the 
prepotential with respect to the Whitham parameters g k . Starting from the expressions 
dHJ) it results 

d 2 T ( x k \ d 2 T ( x k \ 

i?esoo + — duJi , = Res —dQ . (21) 



<-j UU-r I ; l 13 <~U O »T I / 

osidg-fc V k J otodg k V fc 

The right hand side of these formulae express the susy breaking contributions in ()13|) 
as geometrical quantities which can then be read directly as residues. Nevertheless we 
have to remind the interpretation of the mixed second derivatives appearing in (|13p. As 
already mentioned, they should be thought to come from an appropriate higher degree 
system taking gj — > (j > n + 1), with the genus of the curve and (to,Si) kept fixed. 
Using pijl. the residues can be computed directly with gj = 0; therefore, Q to Q g . | fli =o and 
ds dg I g j =0 can t> e properly obtained from the curve of (n+l)-degree which depend only on 
the couplings g k , k — 1, ■ ■ ■ ,n+ 1. We can then extract all the mixed second derivatives, 
included those with respect to gj, using the (n + l)-degree geometry. 
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This simplification is one of the advantages of the embedding of the geometry in 
the Whitham framework. With this approach we compute the holomorphic supersym- 
metry breaking terms in the effective glueball superpotential corresponding to a non- 
super symmetric perturbation of the (n + l)-degree tree-level superpotential (fTTjl. without 
the explicit knowledge of the prepotential T . 



5 Tree— level cubic superpotential 

We consider the simple case of a supersymmetric U(N) gauge theory with tree-level 
superpotential 

7TI Q 

W tree (<S>) = yTr$ 2 + |Tr$ 3 . (22) 

As suggested before we break supersymmetry promoting the coupling constants of the 
tree-level superpotential to spurions ()9I1(JJ) deforming (|2"2"j) as 



Wtreem = m + g2F2 Tr $ 2 + 9 -±^Tr $ 3 + 6 2 V ^Tr & . (23) 
2 3 i>3 J 

The geometry of the M = 1 solution is described by the following complex curve of genus 
one with meromorphic differential 

y 2 = g \ x - ai )\x - a 2 f + fo + hx , (24) 
dS = y dx = \/g 2 (x — ai) 2 (x — (X2) 2 + fo + 2gt x dx . (25) 

with ai = and a 2 = — y the classical roots. 

We do the computation of the supersymmetry breaking parts as a series with small 
width of the cuts and then small values of Sj and to- The approach is the same as in j2]. In 
particular, we have considered the case of classical susy vacua with unbroken gauge group 
and also the case with U(N) — > U(Ni) x [/(A^) gauge symmetry breaking pattern. Using 
fl2"T]) . we compute directly the second mixed derivatives of the prepotential, i.e. the susy 
breaking contributions. The details of the computations are in appendix A. We express 
our results directly in terms of the physical glueball superfields Si {i = 1, . . . , n) using the 
change of variables (jHJ) at the end of the computation. We will write explicitly only the 
novel supersymmetry breaking contributions to the low-energy glueball superpotential 
referring the reader to the literature [3 for the well known M = 1 part. 

In the case U(N) — > U(N), using (jTSj) with superpotential of the form we find 



^Weff = -±W$f{S,m,g) + 



9 2 T, 




T(f ) (8g 2 S 



< (k + l)\ r(|) V rn 



+ 



8 



+ e 2 r 3 



+ e 2 v A 




(k + iy. r(|) 



8g 2 S 



+ 



m 



64o 4 ^ fc! 

3 A:=2 



r(U') 



i>4 J 



5 / # 



T(i(A; + l)) J \m 



8g 2 S 



+ 



i ! W V(9 + ™Z) 2 + f ?-4gSe 



5=o 



2(J - 1)! 



[l + Y) 



dV' 1 ^(g + m0 2 + he-±gse 



(26) 



Y and / are functions of (S 1 , m, g) whose expressions (fo"21 are given in appendix A. 

As a consistency check of our computation and focusing on the spurionic terms T 2 and 
T 3 , we can compare the previous result with the mixed second derivatives of the M = 1 
perturbative prepotential 



s 4z 



r(f) (8g 2 S 



^(* + 2)!r(f + i) 



(27) 



derived for the first time in JE] from the large-N matrix model. We find a complete 
agreement except the linear term (~ S) in the series multiplied by T 2 . 
The appearance of the linear term can be explained in the following way. It is known 
|17| 01 E] that the measure in the matrix model partition function and also the al- 
lowed divergent modes on the complex curve |3j give schematically a contribution like 
(S — S\og(mAl/ S)) where A is a cut-off: this contribution together with the additive 
term (2iriTS) in the effective superpotential gives the Veneziano-Yankielowicz superpo- 
tential [TSJ - The derivatives of this contribution w.r.t. the coupling m give exactly the 
linear term appearing in (|26|) which also agrees with what we have found using perturba- 
tive techniques (see Sec. 6). 

The supersymmetry breaking part coming from the quartic term (and also from the higher 
ones) can be checked by comparison with the M = 1 superpotential computed implicitly 
in for a generic tree-level superpotential. By evaluating the derivative where all the 
coupling constants except (m, g) are set to zero, we find agreement with their computa- 
tion for all the finite order explicitly given by them. 

In the case U(N) — > U(Ni) x U(N 2 ), we consider only r 2 , r 3 , T 4 in as source of 
susy breaking. Then, using (fT5j) we have 



W, 



eff 



Wf f j\S^S 2 ^g) + 



9 



+ 



+ 



+ 



e 2 r, 



(2N 2 - Ni)— + (2iVi - N 2 ) — + 30(JVi - iV 2 



+ 3(5iV 2 - 2^)^-^ + 3(2A 2 - 5iVi)-4S| + O^ 3 



+ 



# 2 r 3 



2N 2 — - 2N±— + 20{N 2 - N 1 )^S l S 2 + 



9 



9 



+ 2 (2 A; 



- hN 2 )-?-Sl + 2(5JVi - 2N 2 )-^-S 2 2 + 0(S 3 ' 

777 777 O 

2N 2 -S 1 + (2Ni + N 2 )-S 2 + -5(2JVi 
^(A 2 - 2^)4 - l(N, - 4iV 2 )4 + 0(S 3 



e 2 v, 



+ 



+ 



3iV 2 )Si,S 2 + 



(2f 



where we show terms up to the quadratic order in S; we give in Appendix A a sketch of 
the computation. 

We can check also this case using the results of |3] for a quartic tree-level superpo- 
tential. Taking the derivatives of their results with respect to the coupling constants and 
then making the appropriate limit (5*3 = and g± — > 0) we get exactly our supersymmetry 
breaking contributions. 

Observe that, also in this case, linear terms appear in the supersymmetry breaking se- 
ries multiplied by T's. These can again be understood as coming from the Veneziano- 
Yankielowicz piece of the effective superpotential. In fact, the scales A< associated with 
each unbroken gauge group sector U(Ni) are functions of the coupling constants as a 
consequence of the threshold matching [20]; by taking derivatives w.r.t. the couplings we 
get exactly those linear contibutions appearing in (|28p. 

Finally we note that, up to the quadratic order in S = Si, we can consistently get our 
first result (J2fij) from the second one (J28j) simply by setting (S 2 — 0, N 2 — 0). 



6 Perturbative arguments 

In this section we exploit the perturbative approach 4J to discuss, from a field theoret- 
ical point of view, our use of the Af = 1 prepotential to study the low-energy glueball 
superpotential in the case with broken susy. We consider only the case of unbroken 
U(N) gauge group and tree-level superpotential for the adjoint chiral superfields given 
by W($) = Ylk=2 ^Tr $ fc where G k = Qk + 2 T k are the spurionic coupling constants. 

We recall that, because of holomorphicity, in the Af = 1 case the effective superpo- 
tential is a function only of the coupling constants gj~ and not of the ~g k |!4[ HJ In our 
case susy is broken by the spurions and holomorphicity in the couplings is not any longer 
a property of the superpotential. 

In a perturbative framework the spurions G k can be thought as ordinary background 
chiral superfields. We can then think susy unbroken and the perturbative computations 
in a superspace approach go using the usual D-algebra |2T] . The effective action will be 
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schematically of the form 



J d 2 9d 2 9 lC{G k , G k , D 2 G k , D 2 G k , ■ ■ ■ , S, S) + J d 2 9 W eff (G k , S) + h.c. , (29) 

where the superpotential W e ff(G k , S) is constrained to be a holomorphic function 7 of G k . 

If we choose the particular supersymmetric configuration in which all the chiral super- 
fields G k are equal to the constants g k , without any dependence on 9 a , then the M = 1 
effective superpotential is W e ff(g k , S) and hence its holomorphicity [Hj. 

If we choose instead the configuration G k = g k + 9 2 T k (T k ^ 0), we break susy and 
furthermore we will have two kind of contributions to the glueball superpotential. 
The first ones come from W e ff(g k + 9 2 T k , S) in (}2*9"j) and are the holomorphic ones we have 
studied in the previous sections. We call them the holomorphic contributions. 
The others are D-terms contributions holomorphic in 5* but not necessarily in the coupling 
constants g k , ~g kl T k and T k which come from particular contributions to K, in ()29)1 and 
which can be written as J d 2 9 integrals contributing to the glueball superpotential 8 . These 
terms in the M = 1 case (r — > 0) are zero. 

Here we adopt a pragmatic attitude and we study only those contributions to the 
glueball superpotential which can be computed using the powerful perturbative techniques 
developed in [I] for the M = 1 case. 

In j3] the perturbative series was generated using only the propagator of the chiral matter 
superfield sector and the antichiral superfield $ was integrated out. This was the central 
point for their simplifications. In order to be able to integrate out 5> as in |3] we must have 
interactions only in terms of the chiral superfield $ and then we consider the following 
UV action 

S($, $) = J d 4 x d 4 9 Tr e- y $e y $ - J d 4 x d 2 9 ^Tr $ 2 - J d 4 x d 2 9 ^Tr l> 2 + 

+ I d 4 xd 2 9 i(# 2 r 2 )Tr$ 2 + I d 4 xd 2 9 jr^-(g k + 9 2 r k )Tr$ k , (30) 

J J 1 Q 



fc=3 



where all the antiholomorphic interactions / d 2 9 [|(6» 2 r 2 )Tr $ 2 + XT=3 f (#fc + 2 I\)Tr 
are neglected. Furthermore, since we are interested in the glueball superpotential it is 
also possible to do the usual simplifications of |U E] finding as the relevant action 9 

J d 4 xd 2 9 j 2=$[n - iW a d a - mm]$ + j , (31) 



We are thinking about the case with masses in the Wilsonian approach for which the nonrenormal- 
ization argument works without IR patologies. 

8 For example, the reader could think about two terms like (with W a — iD 2 (e~ v D a e v ) \21\ ) 
J d 2 6d 2 6G{g,g,Y,T,e 2 )6 2 T S p = fd 2 8GTS p or / <f9d 2 e H (g,g, T,T, 6 2 )Tr[i( e - v D a e v ) W a }TSi = 
Jd 2 9HTS q+1 where G and H arc functions of 5 ,i;,r,r,6> 2 and not 9 2 . 

9 W a = [W a , • • • } is the spinorial gauge field strength adapted to the action, as a graded-commutator, 
on the adjoint representation of the U(N) gauge group. 
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where W™e e m our sus y broken case consists in the second line of (J3UJ). The difference 
with respect to [I] is that the tree-level superpotential is now defined in terms of spurionic 
coupling constants. 

Now, from (pHj) . it is clear that the glueball superpotential we are going to compute 
will be holomorphic in S and in all the coupling constants except, at most, for the mass. 
In particular we observe (we refer to Appendix B for the details) that we can have con- 
tributions only of the following form 

J d 2 6^W eff {G k ,S) + ^J2 B ^k,^k,0)S 1 ^ . (32) 

W e ff(gk + 2 Tk, S) is the holomorphic contribution we have already defined. Instead, the 
second part of (}3*2~|) is a particular subclass of the D-term contributions discussed before 
where B\ are holomorphic in all g k , T k and possibly depend also on 9 2 . 

A careful perturbative analysis of (|3*2*)l shows that all the coefficients £>/ = vanish 
V/ and that the contributions to the glueball superpotential we are computing have the 
following form 

with F = J2 :F °,i(9k + 2 r k )S l . (33) 
i 

We refer the interested reader to Appendix B for the technical details of our perturbative 
computations. 

In (J3*3j) JF 0) z are the planar amplitudes with / index loops of the dual matrix model H] 
where the coupling constants are in this case the spurions G k = Qk + 2 T k - The first term 
in (|3*3*|) is given by a 1-loop diagram with one vertex I^I^Tr^ 2 . This term is associated 
with the 1-loop matter contribution to the Wilsonian beta function for the gauge kinetic 
term which is implicit in the nonperturbative Veneziano-Yankielowicz superpotential. In 
the previous section we have seen that this term is also given by the geometrical methods. 

We conclude that, within our stringent approximations and in the case of unbroken 
U(N), the effective glueball superpotential in the presence of spurions (jHHJ) can still be 
deduced from the Af = 1 holomorphic superpotential supporting the results of the previous 
sections. 

7 Conclusions 

The Dijkgraaf-Vafa conjecture with supersymmetry breaking is the subject of this work. 
We have considered the simple case of U(N) gauge theory with massive adjoint chiral mat- 
ter multiplet with a polynomial tree-level superpotential. We have studied the case where 
supersymmetry is broken in the tree-level superpotential by promoting the coupling con- 
stants to chiral spurions. We have considered their F-components as non-supersymmetric 
small perturbations of the Af = 1 gauge theory and we have discussed how holomorphy 
can still play a role. The non-supersymmetric holomorphic contributions to the effective 



d 2 9 N6 2 T 2 - + N^ 
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low-energy glueball superpotential have been derived with geometrical methods embed- 
ded in the Whitham framework as well as with techniques of superfield formalism with 
spurionic fields. 

Non-holomorphic .D-terms, soft breaking via gaugino mass, low energy vacua are open 
to investigation. This goes beyond the information encoded in the holomorphic matrix 
model that we have used so far. 
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A Solution of the broken superpotential 

In this appendix we show the main tools and details used for the computation in the cubic 
tree-level superpotential of section El 

We have already written the genus one Riemann surface characterising the solution 
(|24 | 12 5 |) . We observe that there is one holomorphic differential ^ defined on this surface. 
This differential can be expanded around the point at infinity in powers of £ = l/x 



, (34) 

6=0 



where R m are functions of to, fo and can be simply computed by power expansion of 
y. The normalized holomorphic differential duj is then 

du = - , (35) 

h y 

where we have introduced the following quantities 10 

h m = (f> . (36) 

y 

The meromorphic differentials dQ k are defined by 11 

cWS xdx ldfodx ddS 
dQ = = 9 + 77^7 > dn k = ^— k = 2, 3 , (37) 

and are completely fixed by the normalization constraints 

dtt = , (h dtt k = k = 2, 3 . (38) 



Then df2 results to be 

in f h\\ dx 

dil = [gx - g— — , (39) 
V ho J y 

where we used the first normalization condition of (f3*K|) which implies |^ = —2g^. 
Collecting these formulas we can express the second derivatives of the prepotential char- 
acterizing the susy breaking terms in the effective superpotential (fT3*j) for the case under 
consideration as follow 

= Rea ° (x dn °) = s ( 9Rk ~ % Rt - 1 ) ' m 

B". OA = %4 , (41) 



ds2dg k \ k J k ho 

10 The a-cycle encircle counterclokwise the second cut accordingly to our conventions 
11 We denote g^ = m and 33 = g of 
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where the Rk are defined in 

In the case of unbroken gauge group (U(N) — > U(N)) the cut associated with the s 2 
variable degenerate to a point with s 2 — > and the only variable is to- The curve (124)1 
can be written as 

y 2 = g 2 (x - x 1 f(x - x 3 )(x - x 4 ) , (42) 

where x 3 , x 4 are the extremal points of the first cut and X\ is the double zero of the curve 
where the second cut degenerates. Following [2] it is useful to introduce the quantities 

A 43 = ±(x 4 -x 3 ) , A = ( ai -a 2 ) = 'f , (43) 

Q = fat + x 3 + 2x 1 ) = (a 1 + a 2 ) = -* , (44) 

/ = |(x 4 + x 3 - 2x x ) = y/A 2 - 2A1 3 , (45) 

x x = 2^ , a = ^ , a = 8aS. (46) 

The above relations can be proved comparing (|24J) and (j42|) . We have also directly intro- 
duced the physically relevant variable S = Si = — y. 

Being interested in finding as in (j4~Uj) we have evaluated hi/ ho in this case 



Then we find 



hi Q-I tA7 , 

ir a =Xl = ~- (47) 



dt dgk k \ ho J k 2k 
where Y is defined as 

Y = i = f I W- (49) 

Then 7 = I(Af 3 ) is a function of A| 3 . Written in terms of X{ the variable to is 



to = —ReSoo(dS) = —Res OQ [g(x — x\)\/ (x — x 3 )(x — x 4 ) dx\ = 

= ^(2xi -x 3 -x 4 )(x 4 -x 3 ) 2 = -|a 43 I . (50) 

From |3DD we find 

<7=(1-Y 2 )Y, (51) 

which gives Y, and then g^f^, as a function of a = 8aS. Solving (J5TJ) and taking the 
appropriate branch we obtain 

i 

Y _ gj (\/^-^-^) 3 l^(8a5) T(i(3t-l)) 

' ^TZ-.V 3 = *h « r<# + D) ■ (52) 
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Once Rk (k = 1,2,3) are found from (|34|) we have the first two softly broken terms of 



d 2 T S 



dt dm m 



1 + 3 



^ (8aS) k r(f ; 



+ r(|; 



(53) 



d 2 T 2S ^ (8aS) fc 



cM 5 5 ^(A; + l)! T 



2 



For the third term of (|2*6*j) the computation is a little more involved. The derivative to be 
computed is 

1)1 T 9 -R, + ^R 3 (l + Y) . (55) 



dt dg 4 4 
The coefficients R3 and -R4 are 

-R4 = 77^7 H 3-5 r , i?3 = — r , (56) 

2# 3 r # 5 £ 4 g 2 



where the unknown function f appears. To compute it we integrate in S the equation 
that can be obtained from the first constraint in ()38|) and from ()47|) 

dT =- 2g h- = ~ 2g -^ = m + gI - (57) 



We then have 



MS] = -2mS - 2 9 [ ™Y[S\dS = Cl - 4™S + f, £ < 8 °f t ^ = 4)) , ( H 



where Ci is a function only of the couplings m and g. Using the other constraints in (|38|) 
that define the derivatives of /o with respect to the couplings (m, g) it can be proven that 
Ci vanishes. Finally, we use the formulas (J551 loTH with ci = and find 

* (8a^ / r(|(3fc - 4)) _ r(j(3fc-l)K 

^0^4 ev w v [ 1 r(ifc) 4 r(i(* + i)) J " 1 j 

We observe that with these ingredients one has formally all the needed quantities to com- 
pute in a closed form, as power series of S, all the susy breaking terms in the effective 
superpotential (|13|) coming from higher order supersymmetry breaking deformation in the 
tree-level superpotential (jlljl - These are functions only of Y(S) and fo(S) as shown in 
P5|). In fact the coefficient are defined as and are functions only of to = — 2S, of 
the couplings, and of fo(S) which is known (|55j). At the end the result is 



The computation in the case of broken gauge group (U(N) — > U(Ni) x U(N 2 )) uses a 
procedure as 0, making the calculation as a power series in the width of the cuts. Our 
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aim is again to compute the susy breaking terms appearing in (|13|) using the formulas 
()4U1 |4"T]) . The main difference with the unbroken gauge group case is that now the curve 
does not degenerate 

y 2 = g 2 (x — Xi)(x — x 2 )(x — x 3 )(x — x 4 ) . (60) 

We introduce quantities analogous as before 12 

A 43 = 5(2:4 - X3) , A 2 i = ~(x 2 - xx) , (61) 
Q = ^{x 4 + x 3 + x 2 + xi) = (at + a 2 ) = —j , (62) 

/ = \{x A + x 3 - x 2 - x x ) = VA 2 - 2A 2 3 - 2A| . (63) 

We don't have anymore the simplification ()47|) and we have to write the integrals ho and 
hi as power series in the widths of the cuts (O(A^J). We then find the inverse expression 
of the widths of the cuts A ab as a functions of (s 2 , t ) and obtain (h , hi) in terms of 
(«2, *o)- 

We have also to evaluate the parameters Rk- They have the form (|54*j) but now f has 
to be understood as a function of two variables and t Q = —2 (Si + S 2 ). Precisely f is a 
function of to an d s 2 which are the independent variables and it is determined through 
the relations 

dto 29 h ' ds 2 fl £$ ho' {M} 

The first equation comes from the normalization condition (|3*K|) while the second one is 
a consequence of the definition (jHJ) of the variable s 2 using the explicit form (|2~K|) of the 
differential dS. 

We then compute directly the second derivatives of the prepotential, the susy breaking 
terms (fTSj). using (|4"01 l4*Tj) . At the end of the computation we change variables (JS} to 
express the superpotential in terms of the physical glueball superfields Sj. What we find 

is (J2HJ). 



B Details on the perturbative approach 

In this appendix we explore the details which give ()33|) from f)3Uj) . We will use and extend 
to spurions the method developed in |3J E2J also reviewing, for reader convenience, their 
basic steps. 

Starting from ()30|). the propagator is the same as in 0] 

— 77? 

< $$ >= tt— . (65) 

□ - mm - iW a d a 

The gauge field strength is considered constant then the bosonic and fermionic integrations 
completely decouple in the computation To compute contributions to the glueball 

12 Note that as concern the classical roots there are no modifications. 
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superpotential, we will use the usual chiral-ring properties of W a |H|3]. For example, we 
will use Tr(W a ) n = with n > 2. 

Using the double line notation, a Riemann surface (oriented for U(N)) with genus g 
is associated to each topologically relevant diagram with L momentum loop and I index 
loop, so that L — I + 2g — 1. The .D-algebra is exactly as in [I]. The only difference is 
that performing the .D-algebra some d a can act on the 9 of the spurions G% = gk + 9 2 Tf s . 
giving new terms. 

It is possible to do some general considerations using the constraints given by the D 
algebra structure, the properties of W a and the geometry of the diagrams in the ampli- 
tudes. 

We fix a diagram with L momentum loops and I = L — 2g + 1 index loops. From the 
W a properties it follows that, for a relevant amplitude for the glueball superpotential, the 
maximal number of allowed W a is 21 otherwise we would have at least one index loop with 
more than three W a . Furthermore, in order to perform the fermionic loop integrations it 
is necessary to have at least 2L d a and then at least 2L W a . The number of W a (#W a ) 
in a non-trivially zero diagram then satisfies the inequality 

2L < #W a <2l = 2L + 2-Ag . (66) 

This implies that g = and the only relevant diagrams to be considered are planar. 
Moreover the relevant contributions to the glueball superpotential have #W Q = 2L, 2L+2. 

We consider first the case jfW a = 2L i.e. #<9 a = 2L. In this case the .D-algebra has 
to be done only inside the fermionic loops and then no derivative acts on the background 
spurions. This is equivalent to say that these contributions are insensible to the 9 de- 
pendence of the Gk- Then for these kind of terms we can reabsorb the quadratic vertex 
!# 2 r2Tr ( 3> 2 into the propagator ()65|) by simply doing the redefinition m — > m + 9 2 T2- 
This is clearly true exept a 1-loop amplitude with one vertex |6 |2 r 2 Tr $ 2 contracted with 
the propagator (|65|) which gives the first term of The resulting contribution to the 
glueball superpotential with #W a = 2L is then given by and, except the linear term 
in S, these term are computed perturbatively using the dual matrix model of the M = 1 
case. 

Now, we consider the other case j^W a = #<9 a = 2L + 2. All index loops are now 
saturated and we have a contribution proportional to S^ L+1 \ 

These contributions are nonholomorphic since they are proportional to m~ 2 . In fact, 
the corresponding diagrams have a multiplicative factor m Pi from the numerator of the 
Pi propagators (|63j). Expanding the propagators (|B3j) at order (W a d a Y 2L+2 \ we have 
Pf = Pi + 2L + 2 bosonic propagators 2+ mm expressed in momentum space. By redefin- 
ing the bosonic loop momentum variables p 2 fnp 2 , from the bosonic Jacobian we are 
left with a contribution m 2L while from the denominator of the bosonic propagators we 
have a term ffi~( Pi + 2L+2 \ Summarizing we have [m Pi ][m 2L ][m~( Pi+2L+2 - ) ] = m~ 2 . There- 
fore, along the calculation we will set fn = 1 and multiply the final result by WT 2 . 
In performing the calculation it is convenient to express the propagator ()65|) in the 
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Schwinger variables 

POO 

I dsi exp [—Si(pi + iWfd a +m)\ 
Jo 

As in the bosonic contribution is given by 
1 1 



(67) 



J boson 



(4tt) 2L (detM(s)) 2 



M ab (s) = ^ SjL ia L ib , pi = ^2L ia k a . (6£ 



We note that M(s) is an L x L matrix and then the denominator of Zb oson (168)1 is a 
homogeneous polynomial of degree 2L in S{. Furthermore, we have, from the fermionic 
integrations of these diagrams, 2L + 2 SjWf terms. Then, at the numerator we have 
a homogeneous polynomial of degree 2L + 2 in Sj. The degree in Sj of the numerator 
results to be greater than the denominator degree. Thus, for the class of diagrams with 
ifW a = 2L + 2 (certainly) there is no cancellation between the bosonic and fermionic 
integrations in contrast with the case jfW a = 2L 

Performing the .D-algebra we realize that there are two distinct possibilities depending 
on the way the two extra d a are distributed on the external spurionic terms. The first 
possibility is when two d a act on one spurionic constant 6 2 T k . This contribution would 
have a multiplicative factor [(1^) n^T^ (S'ku + ^ 2 ^k v )} (V is the number of vertices of the 
considered diagram). 

The second possibility is when the two d a act on two different spurions. In this case we 
have a multiplicative term [(8 a T k )(8 a T k <) Yl^i (dk v + ^fcj]. 

Summarizing the previous considerations, the general structure of the glueball super- 
potential, due to the integration of the matter fields considering only the holomorphic 
part of the interaction vertices as in (jSHJ), is 



<( m 2 r 2 - + n V r Qjl {g k + e'r^is 1 - 1 + -^ V b^, r k , e)s l \ . < 69 > 

m z — ' m A — ' 

i i ) 



The Bi are holomorphic in all g k , and possibly depend also on 6 2 . Bi are analytic in 
all the variables except m. We will show that Bi — V/ justifying (jHSJ)- 

To compute B% we must perform the D-algebra and treat the group theoretical factor. 
As in |1] to simplify the fermionic integrations we can use the fermionic Fourier momentum 
representation. The novelty in the computation is due to the fact that now there is also 
the 9 2 from the spurionic vertices to be Fourier transformed. In particular, we have 



6 i = -^ 2 \6) = - j d 2 7re l 



(70) 



We focus on a planar diagram with L = 1 — 1 bosonic momentum loops with P propagators 
and V vertices. In particular we consider the case with only one spurionic constant 
on which the D-algebra acts nontrivially. The other case is analogue. 
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The fermionic contribution results to be 13 
v P 



J fermion 



v=l i=l 

i=l v=l i=l v=2 

« P-V+2 ki P 

d 2 6 / J] d 2 n a d 2 r 6^ ( 7tj 1 + r) J] [ e -*^*] 

•' 0=1 j\=\ 1=1 

d 2 6 j f[d 2 K a f[[e- s ^>">] . 

J o=l 1=1 



e i(££ =1 *iv) a *° 



(71) 

In the second line 7tj v are the spinorial momuntum connected to the v-th ^-vertex and sat- 
isfy 7Tj v = ±7Ti where the sign is + (— ) if the spinorial momentum is going outside (inside) 
the vertex v. In the last two line we exploit the relations X^J=i 7r 3 -„ =0 {v = 2, • • • , V — 1) 
with which we define the remaining I spinorial variables K a from the independent 7Tj. Fur- 
thermore, in the last line we have used the fact that we are searching for a contribution 

to Bi which has = 21 = 2L + 2. Then, in the expansion of e~~ j v =injv a we can 

keep only 1, the term independent of 9 a . This is equivalent to say that only the term in 
which 9 2 of the spurion is killed by two d a contributes to T&\. 

At the end of the above manipulations we remain with / fermionic integrations over the 
indipendent variables n a and the 7Tj are linear combinations of them 

i 

T^ia = ^ Li a K aa . (72) 
0=1 

As in 0] we can implement the requirement of having two insertions of W a for each index 
loop introducing 21 auxiliary grassmanian variables adapted to the action on the 
adjoint representation with 

i 

771=1 

The matrix K is defined so that for each oriented z-th propagator the m— th index loop can 
coincide and be parallel giving K im = 1; or coincide and be anti-parallel giving K im — — 1; 
or not coincide giving K im = 0. 

13 The index j v , depending on v = {1, • • • , V}, runs from 1 to k v which is the degree of the interaction 
vertex v as: Tr $ fc " . 
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Summarazing we find from the fermionic integration 




aa 



) 



(16it 2 S) 1 / U d 2 K a d 2 W rn exp -Y,WlN{s) ma K, 




^(47r) 2 '(detAT( s )) 2 , 



(74) 



with 



N(s) ma = J2 s i K LL 



(75) 



The relevant fact is that, for our class of diagrams which has an S 2 topology, the matrix 
K has a nontrivial kernel. In fact, for example, the vector b m , whose components are all 
equal to one, belong to the kernel of K im 141 . This is simply due to the fact that in the 
case we are studying all momentum propagator lines have exact two index loop passing 
through them with opposite orientation; then, Vi there will be only one Ki m > = 1 and 
one K im ii = — 1 (m! ^ m") and ^2 m K im b m = 1 — 1 = 0. It follows that also the matrix 
Wi s )}am = 12i s i^iK im has a nontrivial kernel indipendently of the explicit form of L 
which we have not analyzed in detail. Then det (iV(s)) = 0. This imply that B\ = VZ as 
claimed before. 

References 

[1] R. Dijkgraaf and C. Vafa, "Matrix models, topological strings, and supersymmetric 
gauge theories," Nucl. Phys. B 644 (2002) 3 [arXiv:hep-th/0206255| . 
R. Dijkgraaf and C. Vafa, "On geometry and matrix models," Nucl. Phys. B 644 



(2002) 21 |arXiv:hep-th/0207106| . 

R. Dijkgraaf and C. Vafa, "A perturbative window into non-perturbative physics", 
[arXiv:hep-th/0208048| . 

[2] F. Cachazo, K. A. Intriligator and C. Vafa, "A large N duality via a geometric 
transition", Nucl. Phys. B603 (2001) 3, |arXiv:hep-th/0103067| . 
F. Cachazo, B. Fiol, K. A. Intriligator, S. Katz and C. Vafa, "A geometric unification 
of dualities," Nucl. Phys. B 628, 3 (2002), |arXiv:hep-th/0110028| . 

[3] F. Cachazo and C. Vafa, "N = 1 and N = 2 geometry from fluxes," 
|arXiv:hep-th/0206017| 



14 Our susy broken case is similar to the situation which appears in the study of the perturbative 
reduction to matrix models for the case of N — 1 supersymmetry and SU/SO/Sp(N) gauge groups 
developed in [55] extending 0]. 



21 



[4] R. Dijkgraaf, M. T. Grisaru, C. S. Lam, C. Vafa and D. Zanon, "Perturba- 
tive Computation of Glueball Superpotentials" , Phys.Lett. B573 138-146 (2003), 
|arXiv:hep-th/02110T7| . 

[5] F. Cachazo, M. R. Douglas, N. Seiberg, and E. Witten, "Chiral rings and anomalies 
in supersymmetric gauge theory", JHEP 12 (2002) 071, |arXiv:hep-th/0211170 . 

[6] L. Girardello and M. T. Grisaru, "Soft Breaking Of Supersymmetry," Nucl. Phys. B 
194, 65 (1982). 

[7] A. Lawrence, and J. McGreevy, "Local string models of soft supersymmetry break- 
ing", JHEP 0406 (2004), 007, |arXiv:hep-th/0401034| . 

A. Lawrence and J. McGreevy, "Remarks on branes, fluxes, and soft SUSY breaking," 
arXiv:hep-th/0401233| 

[8] A. Gorsky, I. Krichever, A. Marshakov, A. Mironov and A. Morozov, "Inte- 
grability and Seiberg- Witten exact solution," Phys. Lett. B 355 (1995) 466 
|arXiv:hep-th/9505035| . 

T. Nakatsu and K. Takasaki, "Whitham-Toda hierarchy and N = 2 supersymmetric 
Yang-Mills theory," Mod. Phys. Lett. A 11 (1996) 157 |arXiv:hep-th/9509162| . 
H. Itoyama and A. Morozov, "Prepotential and the Seiberg- Witten Theory," Nucl. 
Phys. B 491 (1997) 529 |arXiv:hep-th/9512161| . 

[9] L. Chekhov and A. Mironov, "Matrix models vs. Seiberg- Witten/ Whitham theories", 
Phys. Lett. B552 (2003) 293, |arXiv:hep-th/ 0209085|. 

L. Chekhov, A. Marsha kov, A. Mironov and D . Vasiliev, "DV and WDVV" , Phys. 
Lett. B562 (2003) 323, |arXiv:hep-th/0301071| . 

H. Itoyama and A. Morozov, "The Dijkgraaf- Vafa prepotential in the context of 
general Seiberg- Witten theory," Nucl. Phys. B 657 (2003) 53 |arXiv:hep-th/0211245| . 

[10] L. Alvarez- Gaume, J. Distler, C. Kounnas and M. Marino, "Large Softly broken 
N = 2 QCD," Int. J. Mod. Phys. A 11 (1996) 4745 |arXiv:hep-th/9604004| . 
J. D. Edelstein, M. Marino and J. Mas, "Whitham hierarchies, instanton corrections 
and soft supersymmetry breaking in N = 2 SU(N) super Yang-Mills theory," Nucl. 
Phys. B 541 (1999) 671 larXiv:hep-th/9805172| . 

[11] N. Seiberg and E. Witten, "Electric - magnetic duality, monopole condensation, and 
confinement in N=2 supersymmetric Yang-Mills theory," Nucl. Phys. B 426 (1994) 
19 [Erratum-ibid. B 430 (1994) 485] |arXiv:hep-th/9407087| . 

[12] F. Cachazo, N. Seiberg and E. Witten, "Chiral Rings and Phases of Supersymmetric 
Gauge Theories", JHEP 0304 (2003) 018, |arXiv:hep-th/0303207| . 

[13] D. Kutasov, A. Schwimmer and N. Seiberg, "Chiral Rings, Singularity Theory and 
Electric-Magnetic Duality," Nucl. Phys. B 459 (1996) 455 |arX iv:hep-th/9510222| . 



22 



[14] N. Seiberg, "Naturalness Versus Non-Renormalization Theorems", Phys. Lett. B 
318 (1993) 469, |arXiv:hep-ph/9309335| . 



[15] H. Fuji and Y. Ookouchi, "Comments on effective superpotentials via matrix models," 
JHEP 0212, 067 (2002), |arXiv:hep-th/0210148| . 

[16] E. Brezin, C. Itzykson, G. Parisi, and J. B. Zuber, "Planar Diagrams", Commun. 
math. Phys. 59 35-51 (1978). 

[17] H. Ooguri and C. Vafa, "Worldsheet derivation of a large N duality," Nucl. Phys. B 
641 (2002) 3 |arXiv:hep-th/0205297| . 

[18] G. Veneziano and S. Yankielowicz, "An Effective Lagrangian For The Pure N=l 
Supersymmetric Yang-Mills Theory," Phys. Lett. B 113 (1982) 231. 

[19] F. Ferrari, "On exact superpotentials in confining vacua," Nucl. Phys. B 648 (2003) 
161 larXiv:hep-th/0210135| . 

F. Elmetti, A. Santambrogio and D. Zanon, "On N = 1 exact superpotentials from 
U(N) matrix models," |arXiv:hep-th/0506087| . 

[20] K. A. Intriligator, R. G. Leigh and N. Seiberg, "Exact superpotentials in four- 
dimensions," Phys. Rev. D 50, 1092 (1994), |arXiv:hep-th/9403198| . 
K. A. Intriligator and N. Seiberg, "Lectures on supersymmetric gauge theo- 
ries and electric-magnetic duality," Nucl. Phys. Proc. Suppl. 45BC, 1 (1996), 
larXiv:hep-th/9509066| . 

[21] S. J. Gates Jr., M. T. Grisaru, M. Rocek and W. Siegel, Superspace, or One Thou- 
sand and One Lessons in Supersymmetry Benjamin/Cummings, Reading (1983) 
|arXiv:hep-th/0108200| . 

[22] P. Kraus, M. Shigemo ri, "On the Mat ter of the Dijkgraaf-Vafa Conjecture", JHEP 
0304 (2003) 052, |arXiv:hep-th/0303104| . 



23 



